Event-by-event fluctuations of the kaon to pion ratio in central Pb+Pb 

coflisions at 158 GeV per Nucleon 
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We present the first measurement of fluctuations from 
event to event in the production of strange particles in colli- 
sions of heavy nuclei. The ratio of charged kaons to charged 
pions is determined for individual central Pb+Pb collisions. 
After accounting for the fluctuations due to detector resolu- 
tion and finite number statistics we derive an upper limit on 
genuine non-statistical fluctuations, perhaps related to a first 
or second order QCD phase transition. Such fluctuations are 
shown to be very small. 

PACS numbers: 25.75.-q 

Quantum Chromodynamics predicts that at suffi- 
ciently high energy density strongly interacting matter 
will undergo a phase transition from hadronic matter 
to a deconfined state of quarks and gluons, the quark 



gluon plasma (QGP) [Q. To create and study this state 
of matter in the laboratory collisions of heavy ions are 
studied at the CERN SPS which provides lead (^o^Pb) 
beams of 158 GeV per nucleon. Recent data suggest 
that conditions consistent with the creation of a QCD 
phase transition are indeed reached in central Pb+Pb 
collisions H] . Whereas the position of the phase transfor- 
mation in temperature, energy density and baryon den- 
sity may thus be located we are lacking information as 
to the nature and order of that transition. These might 
be reflected by the presence or absence of fluctuations 
that are characteristic for a first or second order phase 
transition. An early theoretical investigation by Kapusta 
and Mekjian [|j suggested such fluctuations in the kaon- 



to-pion total-yield ratio, due to supercooling-reheating 
fluctuations produced by a predicted large enthalpy dif- 
ference in the two phases. The K/i: ratio was shown to 
fluctuate by about 10%, over the domain of conceivable 
hadronization temperatures, i.e. 140 <T<200 MeV Q. 
This prediction would be experimentally testable if the 
K/i: ratio or related quantities could be quantifled for 
individual central collision events H. 

In a single central Pb-t-Pb collision at 158 GeV per 
nucleon about 2400 hadrons are created |^, permitting 
a statistically signiflcant determination of momentum 
space distributions and particle ratios on an cvcnt-by- 
event basis. Using the NA49 large acceptance hadron 
spectrometer, which detects about 70% of all charged 
particles, we are able to study event-by-event fluctuations 
of hadronic observables. In a previous publication ||7| we 
have presented results concerning the fluctuations of the 
eventwise mean transverse momentum. We showed that 
these fluctuations are very small in central Pb-|-Pb colli- 
sions, and that the data are consistent with a hadronic 
gas in thermal equilibrium [pj . 

In this letter fluctuations in the eventwise ratio of the 
number of charged kaons to the number of charged pions 
([K^ +K~]/['K^ +TT^]) are investigated, yielding informa- 
tion on fluctuations in hadrochemical composition P] and 
on the strangeness-to-entropy ratio H. We shall show 
that, similar to the (jpt) study, the event- by-cvcnt fluc- 
tuations of the K /t: ratio exhibit no significant large-scale 
fluctuation signal. We discuss fluctuations predicted in 
a thermodynamical resonance-gas model, which was pre- 
viously used in a study of transverse momentum fluc- 
tuations pj, and those expected in various microscopic 
descriptions of the collision evolution. 

A detailed description of the NA49 experiment can be 
found in |10|. We used a data set of central Pb-|-Pb col- 
lisions that were selected by a trigger on the energy de- 
posited in the NA49 forward calorimeter. The trigger ac- 
cepted only the 5% most central events, corresponding to 
an impact parameter range h < 3.5 fm. The event vertex 
was reconstructed using information from beam position 
detectors and a flt to the measured particle trajectories. 
Only events uniquely reconstructed at the known target 
position were used. In this analysis particles were se- 
lected that had at least 30 measured points in one of the 
two Main Time Projection Chambers (MTPC) outside 
the magnectic field. 

A cut on the extrapolated impact parameter of a par- 
ticle track at the primary vertex was used to reduce 
the contribution of non- vertex particles originating from 
weak decays and secondary interactions. We estimate 
that about 60 % of such particles are rejected by the ver- 
tex cuts. No further acceptance cuts were made, thus 
maximising the statistical significance of the event-by- 
event particle ratio measurement. 

The particle identification (PID) in this analysis was 
based on the measured specific energy loss (dE/dx) in the 
MTPC using a truncated mean algorithm. Details of the 
dE/dx measurement can be found in |10| and in partic- 



ular in [ni|. As the event-by-event particle identification 
depends crucially on the stability of the dE/dx measure- 
ment with respect to time, event multiplicity and possi- 
ble backgrounds, only the energy loss of the track in the 
MTPC was used in this analysis, making use of correction 
algorithms developed for this detector |11|. In particu- 
lar we observed a significant multiplicity dependence due 
to the high charge load on the TPC readout chambers 
in central Pb-|-Pb events. The procedure we employed 
to correct these effects are based on detailed measure- 
ments of the electronics response using laser tracks. They 
perform a channcl-by-channel iterative correction of the 
raw TPC charge measurements taking into account the 
charge history of sets of neighbouring channels, which 
are coupled via crosstalk effects through the sense wires 
of each TPC readout chamber. These corrections im- 
proved the average dE/dx resolution by about 30% from 
f^ dE/dx /{dE/dx) = 5-6% to the final value of 3.9% for 
central Pb-|-Pb collisions. More importantly, the multi- 
plicity dependence of the dE/dx measurement was re- 
duced by more than 90%, leaving a change of less than 
0.3% in the multiplicity range used in this analysis. 
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FIG. 1. Phase-space distribution of accepted charged par- 
ticles pions (above) and kaons (lower plot). No acceptance 
correction was applied to the kaon to pion ratios reported 
here. 

Within our acceptance window for kaons and pions 
(shown in Fig. 1) the dE/dx resolution translates into 
an average separation of pions from kaons of about 2.1 
sigma and of kaons from protons of about 1.8 sigma. Ob- 
viously, given the accepted multiplicity of about 40 kaons 
per event, it becomes crucial to make maximum use of 
the available information when constructing an estimator 
for the event-by-event kaon-to-pion ratio. In particular 
a simple counting of particles is not possible under these 
conditions. In our analysis we use a PID method closely 



related to that proposed in [Q. 

Combining the information from all events, the mo- 
mentum distributions normalized to unity, Fm{p), were 
determined for each particle species (m = kaons, pions, 
protons, electrons). We also evaluated the normalized 
probability density functions for the truncated mean en- 
ergy loss, fm{P', dE/dx), as a function of particle momen- 
tum for each species. The relative yield of different parti- 
cle species is characterized by parameters Qrm such that 
Xlm ©m = 1- These parameters, with the ratio Qr/Qti 
giving the K/ir ratio of interest here, were determined for 
every single event by maximizing the likelihood function 



L 



HEQ' 



Fmipi) fm{Pi; idE/dx)i)] 



(1) 



using directly the observed momentum pi, and specific 
energy loss {dE/dx)i, for each particle i in the event. 
The use of fixed momentum distributions Fm (p) , imply- 
ing small changes in the particle momentum distribu- 
tions from event to event, is justified by the result of our 
analysis of event-by-event fluctuations in transverse mo- 
mentum m which were found to be in the range of 1% 
or smaller. The resulting distribution of event-by-event 
K/tt ratios is shown in Fig. 2 (points). The shape of 
the distribution can be understood as the result of three 
main contributions: 

Firstly, due to the finite number of particles produced 
and observed per event, the ratio of particle multiplicities 
measured event-by-event will exhibit statistical fluctua- 
tions with a width dictated by the individual particle 
multiplicities. 

Secondly, due to non-ideal particle identification these 
pure number fiuctuations will be smeared by the exper- 
imental dE/dx resolution and the event-by-event fitting 
procedure outlined above. 

Finally, superimposed on the statistical and experi- 
mental fiuctuations we expect to observe any true non- 
statistical fiuctuations. The characterization of the 
strength of these non-statistical fluctuations is the goal 
of this analysis. 

In order to isolate the strength of non-statistical fluctu- 
ations in the event-by-event K/tt distribution, we need to 
establish a reference providing us with the size of the con- 
tributions from finite number statistics and experimen- 
tal resolution, but not containing any further correlations 
between particles. This was achieved using a mixed-event 
reference sample. The mixed events were constructed by 
combining particles randomly selected from different real 
events, reproducing the multiplicity distribution of the 
real events. Only one track of each event was used in 
the mixing to exclude any residual correlations. Within 
our central event sample no further selection on impact 
parameter was made. By construction, the mixed events 
on average have the same kaon to pion ratio as the real 
events, but no internal correlations. Due to the con- 
straint on the overall multiplicity distribution they give 
an accurate estimate of finite-number fiuctuations in the 



kaon and pion multiplicities. Equally importantly, each 
track in the mixed events represents an actual dE/dx 
measurement. 
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FIG. 2. Distribution of the event- by- event kaon to pion 
ratio estimated using a maximum likelihood method (points) . 
As a reference, the same procedure was applied to a mixed 
event sample (histogram). 

The mixed events therefore automatically include all 
the effects of detector resolution. They are subjected 
to the same maximum likelihood fit procedure as the 
real events, allowing a direct comparison with the data. 
Any deviation of the distribution for the data from the 
purely statistical mixed-event distribution would indicate 
the presence of non-statistical fluctuations in the produc- 
tion of kaons and pions. The result of this comparison 
is shown in Fig. 2, showing the distribution of the K/tt 
ratio for data (points) and for mixed events (histogram) . 
The small difference between data and statistical refer- 
ence immediately illustrates that any correlations or anti- 
correlations in the final multiplicities are small. In par- 
ticular no significant number of events with unusually 
small or large strangeness content is present in our data 
sample. 

The observed relative width of the data distribution is 
c^data = 23.27%. The width of the mixed event distribu- 
tion is (Jmixed = 23.1%. A morc detailed analysis of the 
latter ||ll| shows that it is composed of contributions from 
finite number statistics of cr = 15.9% and the dE/dx res- 
olution plus fitting procedure of cr = 16.7% which add in 
quadrature to the observed width. Note that the average 
eventwise K/tt ratio shown in Fig. 2 does not necessarily 
correspond to the true ratio in Att acceptance j^], since 
no corrections for efficiency and y,PT acceptance were 
included in the analysis presented. 

To quantify the remaining difference between data and 
mixed events we define the strength of non-statistical 
fiuctuations as 



^non — stat 



data mixed 



(2) 



In general, processes leading to a correlated production 



of one or the other particle species or to a correlation in 
their multiplicities would result in (Jnon-stat > 0. 
For our data set we obtain 



2.8% ± 0.5%. 



(3) 



The fluctuations observed in the data are very small 
compared to the scale given by the observed two-fold 
enhancement of strangeness in nucleus-nucleus collisions 
relative to nucleon-nucleon collisions |]6| . To a very good 
accuracy we can conclude that the mechanism responsi- 
ble for the enhancement of strangeness is therefore active 
in each central Pb+Pb collision. 

The observed value of Onon-stat also allows us to es- 
tablish a limit on event-by-event fluctuations as a func- 
tion of the frequency of occurrence of these fluctuations. 
In establishing this limit at the 90% confidence level we 
took into account our estimated systematical uncertain- 
ties in determining the event-by-event K/n ratio. As 
mentioned previously, all systematic uncertainties con- 
sidered tend to increase the observed value. According 
to our measurement fluctuations of a relative amplitude 
of (Jnon-stat > 4.0% Can bc ruled out in case all events of 
the sample exhibit the fluctuation pattern, whereas fluc- 
tuations occurring in 5% of the events can only be ruled 
out for Unon-stat > 15.0% (stiU Small compared to the 
strangeness enhancement). 

Finally, we would like to compare our measurement to 
the FRITIOF model of nucleus-nucleus collisions. This 
microscopic model quantitatively reproduces the corre- 
lations observed in nucleon-nucleon collisions and pre- 
dicts a value of cr™"^/ = 9% for central Pb+Pb col- 
lisions. Clearly the fluctuations in the data are much 
smaller. The disappearance of fluctuations when going 
to central nucleus-nucleus collisions can be interpreted as 
evidence for statistical particle production, as opposed 
to the highly correlated processes observed in nucleon- 
nucleon collisions. 

This observation lends further support to the interpre- 
tation of particle ratios in nucleus-nucleus collisions us- 
ing statistical hadronization models [g[, with the added 
information that the statistical distribution is actually 
realized in every individual event and not just on an en- 
semble basis. 

The observed fluctuations can also be compared to 
the results of an equilibrium hadron gas calculation M] 
which includes correlations induced by many-body de- 
cays of resonances into pions and kaons. This calculation 
predicts non-statistical fluctuations (Jnon-stat ~ 2%, in 
good agreement with the data. 

In conclusion we have presented the first event-by- 
event measurement of particle ratios in ultra-relativistic 
heavy ion collisions. Using particle identification by 
dE/dx in the NA49 TPCs the fluctuations in strangeness 
production were studied using the ratio of charged kaons 
to pions. No non-statistical fluctuations are observed and 
we deduce an upper limit of (Jnon-stat < 4.0% for fluctu- 
ations occurring in every event at the Sa-level. The fluc- 
tuations are therefore very small relative to the two-fold 



strangeness enhancement, indicating that the dynamical 
evolution of individual events proceeds in a very simi- 
lar fashion. As in our previous study of transverse mo- 
mentum fluctuations in Pb-|-Pb collisions W, we find no 
evidence of fluctuations that might indicate a strong flrst- 
order phase transition or freeze-out near a proposed QCD 
critical end point. 

On the contrary, the fluctuations observed in Pb+Pb 
collisions are significantly smaller than those expected 
for an indepedent superposition of nucleon-nucleon colli- 
sions. This supports the interpretation of fiavor ratios in 
terms of statistical hadronization models, combined with 
a smooth transition from a possible partonic state to the 
final state hadronic particle composition. In fact, the 
minimal fluctuations expected due to resonance produc- 
tion of the final state hadrons seem to completely exhaust 
the observed strength of non-statistical fluctuations. 
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